A New Expansion of the Heisenberg Equation of Motion 
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We derive a new expansion of the Heisenberg equation of motion based on the projection operator 
method proposed by Shibata, Hashitsume and Shingii. In their projection operator method, a certain 
restriction is imposed on the initial state. As a result, one cannot prepare arbitrary initial states, 
for example a coherent state, to calculate the time development of quantum systems. In this paper, 
we generalize the projection operator method by relaxing this restriction. We explain our method 
in the case of a Hamiltonian both with and without explicit time dependence. Furthermore, we 
apply it to an exactly solvable model called the damped harmonic oscillator model and confirm the 
validity of our method. 



I. INTRODUCTION 
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o 
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^jf^, It is important to calculate the time evolution in quantum field theory in various fields of physics. However, there 
^ is no established method to carry out such calculations. In this work, we study a new expansion of the Heisenberg 
\ equation of motion based on the projection operator method. The projection operator method, which we use in this 

paper, was first proposed by Shibata, Hashitsume and Shingii. ||l[] 
■ The projection operator method has several important characteristics, (i) It can be applied both in Schrodinger 
, and Heisenberg pictures, (ii) There is an arbitrariness in the selection of projection operators, (iii) Equations with 
I ' and without time convolution integrals are treated systematically, (iv) One can use it to calculate in the case of an 
^ , explicitly time-dependent Hamiltonian. The method has been successfully applied to problems in quantum mechanics. 
fSl dl The Mori and Nakajima-Zwanzig methods, which are well-known methods using the projection operator, are 
00 included in this projection operator method. We believe that the projection operator method is more convenient 
> than other common formalisms, e.g., the closed time-path formalism and the Feynman- Vernon influence functional 
' technique. |^ In such formalisms, further approximation is needed to derive an equation of motion without a time 
convolution integral, for example, as performed in the physics of hadrons and the early universe, j^] However, in the 
projection operator method, there is no such problem, because of the above characteristic (iii). Furthermore, since 
this method consists of canonical formalism, it is useful to discuss the meanings of time evolution in quantum field 
theory by comparing with usual scattering theory. 
I ' In a previous paper, we discussed the renormalization of ultraviolet divergences in a certain model, using an equation 
of motion without a time convolution integral derived with the projection operator method, j^] However, there is still 
the problem that we cannot prepare arbitrary initial states, because a restriction is imposed on the initial state to 
carry out the systematic expansion of the Heisenberg equation of motion. In this work, we study a new expansion of 
the Heisenberg equation of motion based on the projection operator method without imposing the restriction. For 
this reason, we can prepare more general initial states, such as coherent states, which may be important to describe 
\2{ ' phase transitions. 

5^ ] This paper is organized as follows. In §2, we explain the formalism in the case that Hamiltonian has no explicit 
time dependence. In §3, the case in which Hamiltonian has explicit time dependence is discussed. Such a case is also 
important for studying physical phenomena, for example, pair creation in a strong external field. In §4, we apply our 
formalism to the damped harmonic oscillator model, which is exactly solvable, to confirm the validity of our method. 
Conclusions are given in §5. 



II. THE HEISENBERG EQUATION WITH TIME-INDEPENDENT HAMILTONIAN 

Our strategy is to study quantum field theory as an initial value problem. In this section, we derive our equation 
in the case that the Hamiltonian has no explicit time dependence. Our starting point is the Heisenberg equation of 
motion, 

±0{t)=^[H,Oit)] (1) 
= iLO{t) (2) 
— > 0{t) ^ e'^(*-*")0(io), (3) 
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where L is the Liouville operator and is the time at which we prepare an initial state. The Heisenberg equation 
contains complete information of the time evolution of the operator, but in general, it is difhcult to solve exactly when 
there are interactions. Therefore it is necessary to make some approximations. For this purpose, we introduce generic 
projection operators P and Q that have the following general properties: 



P^ = P, 
g = 1 - P, 
PQ = QP = 0. 

By using these projection operators, we can carry out coarse-grainings in the time development. From Eq. 
can see that the time dependence of the operators is determined by e*^(*~*o). This yields the equation 



(4) 
(5) 
(6) 



dt 



J,L{t-to) _ QiL(t-ta)^-j^ 



,ii(t-to) 



iP + Q)lL. 



From this equation, we can derive the following two equations: 

d 



dt 
d 

—( 

dt 



iL{t-to)Q _ giL(t-io) 



PiLQ + e'^'^*-^°'>QiLQ. 



Equation (0) can be solved for e'^^' *°^Q- We obtain 



zL{t-to)Q ^ Q^tLQ{t-to) ^ I ^gg»L(s-to)pj^Qg»LQ(t-s) 

= ge^-'^'5(*-*o) 



Jto 



S(t,to) 



-e*^(*-*°)pi](i,io) 



l-E(t,io)' 



where 



J:{t,to)= I dse-'^^'-'^PiLQe' 

Jto 

Substituting Eq. ( pi] ) into Eq. and operating with 0(to) from the right, we obtain 



(7) 

(8) 
(9) 

(10) 

(11) 
(12) 



dt 



0{t) 



_ JLit-to) 



PiLOito) + e 



iL{t-to) 



P^{t,to) 



l-E(t,io) 



tLO{to) 



+Qe 



iLQ{t-to)_ 



lLO{to) 



(13) 



This equation is equivalent to the Heisenberg equation, and it has no time convolution integral. We demonstrate 
this point at the end of this section. Therefore, we call this the time-convolutionless (TCL) equation. It is the same 
equation as that used in the previous projection operator method. Q 

Now we consider the situation that the total system can be divided into two parts: the system and the environment. 
We want to know the detailed behavior of the degrees of freedom of the system. Therefore, we carry out coarse- 
grainings to treat the environment degrees of freedom, as in the Feynman- Vernon influence functional technique |^ 
and the projected effective Hamiltonian. |^ Our next task is to specify the projection operator for this purpose. In this 
case, the total Hamiltonian can be divided into three parts, the system (S), the environment (E) and the interaction 
(I) between the system and the environment: 



H = Hg + He + Hj. 



(14) 
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The self- interaction of the system and/or the environment can also be included in Hj. Furthermore, we assume that 
the initial density matrix p is given by the direct product of the system density matrix ps and the environment density 
matrix pE, all given at the initial time t^: 

p = ps®PE- (15) 

We then define the projection operator as 

PO = 1:ye[peO] = {0)e (16) 

for any operator O. With this projection, we can replace an operator acting on the environment with a c-number. 
From the nature of the projection operator, we obtain the following relations: 

PLs = LsP, QLs = LsQ, (17) 
LeP = 0, LeQ = Le, (18) 

where LaO = [Ha,0] for a = S,E,I. Using these properties of the projection operator, we have the following 
relations: 



QLoQ = QLo, (19) 
Qe^^°*Q = Qe'^^"^*Q = Qe'^°\ (20) 

where Lq — Ls + Le- 

In the previous projection operator method, the additional relation PLe = is assumed to derive a systematic 
expansion. This relation is satisfied only in restricted cases, for example, when pE is an eigenstate of He or a 
mixed state which includes diagonal components in the basis of eigenstates of He- Therefore, this relation is not 
satisfied for general pE- Now, we take a coherent state for as an example: 

PS = e-|"l'e"'''|0)(0|e"*" = |a)(«|. (21) 
The state \a) is an eigenstate of the annihilation operator a: 

a|a}=a|a). (22) 

The quantity PLeO is calculated as 

PLeO - TteIpeLeO] 

= {a\[HEO - OHE]\a) 

^ Ea*{a\aO\a) - Ea(a\Oa''\a), (23) 

where we assume He = Ea^a. The r.h.s. of the above equation does not become zero for an arbitrary operator O. 
Therefore, we can see that the relation PLe = is not satisfied. From these facts, we can conclude that the expansion 
used in the previous projection operator method is not applicable to general initial states. We will derive a systematic 
expansion of the Heisenberg equation of motion without imposing the relation PLe — 0. With this, we can calculate 
the time evolution by using more general initial states, but with the assumption Eq. (p^). 

With the above properties of the projection operators (19h and (EQ), the function I](i,to) can be expressed as 



S(i,io)= / dse''^^^-''^PiLQ 

J to 



^iLQ[t~s) 



— Q — g^-iL(t-to) Q^iLQ(t-to) 

= Q - e-''^°'^^-*°^C{t,ta)QV{t,ta)e'^°^*-*°'> 

= Q - e-'-^n(*-*°)Qe*-^o(*"*o) - e-'^°'-*-^°\C{t, to) - l)Qe'^°'^^-*°^ 
-e-'^°^*-'°'>Q{V{t,tn) - l)e'^o{t~to) 

_g-«Lo(t-to)(c(f,to) _ l)Q{V{t,to) ~ l)e'^oit-to) _ (24) 



Here, the functions C{t,tQ) and V{t,to) are expressed as 
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rt /-ti rt,^-l 



^ PL PL I /'Cti — 1 

xLi{tn - to), 



(25) 



'D{t,to) = (,iQLQ{t-to) ^-iQLoQit-to) 



1 + dti dt2--- dtnQLf{tn-ta)QLf{tn-i-ta) 

^ — 1 <^ tn J tr\ J tf] 



f to ^ to J to 



(26) 



where 



Li{t - to) = e*-^«(*-*«)L/e-*^°(*-*°), (27) 
L^(t - tn) = e*-^«(*-*«)i/Qe-*^°(*-*°). (28) 

The operator C{t,to) [P(t,io)] is a time [an anti-time] ordered function of Liouville operators. These expressions are 
derived in Appendix A. Then, we have 

PI](t, to)- = P^it, to) ^ 



l-E(t,io) ' "'l-QS(t,to) 

Pe-*-^«(*-*«)C(t, io)Qe'-^«(*-*") 
1 



1 + e-'-f'«(*-*«)(C(i, to) - l)Qe^^o{t-to) ' 

(29) 



Note that the second hnc on the r.h.s. does not include V{t,to)- The detailed derivation of Eq. ( p9[) appears in 
Appendix B. Substituting Eq. ( p9| ) into Eq. (13), we obtain the final expression of the Heisenberg equation: 

^O(i) ^ e:^^'~'°^ PiLOito) 

-e^Lit-to)p^-^Loit-to)c(^^^^^^Q \ ^^ e'^"^*-'°\LO{to) 

i + (L{t, to) — ijC^ 

^Qg.LQ(t-to)__^__^^Q(^^)_ (30) 

When we expand PJ^{t,to)/{l — I]{t,to)) up to first order in the interaction Hj, we have 

—0(t) = e'^'-*'^°^ Pe-'^«^^-*°'^ Pe'^°'-*-*°hLO(to) 
dt 

_^^rL{t-to) p^-^Lo{t-to) p f dse'-^°("-*«)iL/e-'^°("-*")Qe*-^''(*-*")iLO(to) 



+Qe'^Q(*-*«)^-^^^LO(to). (31) 



Here, we do not expand the third term on the r.h.s. of Eq. (31). This term becomes zero when we take the expectation 
value, and therefore we do not expand it. This expanded equation is used in §4. 

It can be seen that Eq. (|3|) does not contain a time convolution integral, because of the form of the full time- 
evolution operator, 6'^^*"*"^ which operates from the left in the second term on the r.h.s. of the equation. If this 
did contain a time convolution integral, the form of the full time-evolution operator must be 6*^*^*^"^ where s is an 
integral variable. Such a time-convolution (TC) equation is discussed in Appendix C. In other formulations, e.g., the 
closed time-path formalism and the Feynman- Vernon influence functional technique, the derived equation of motion 
has a time convolution integral in general. Such a time convolution term is called a "memory term" . The derived 
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equation of motion is often solved using the Markov approximation. |^ However, in our improved projection operator 
method (and the previous projection operator method), the equation without a time convolution integral, that is, the 
TCL equation is automatically obtained and it is not necessary to make the Markov approximation. 

We would like to make some remarks regarding the operator O{to). First, the choice of the operator O{to) is 
not restricted to an operator of the system. We can choose an environment operator or a product of system and 
environment operators as the operator O(to)- This is different from the case of the previous projection operator 
method and the Uchiyama-Shibata (U-S) projection operator method, which was proposed recently, H because those 
method use the condition PO{to) = O(to)j which is satisfied only when O{to) is a system operator. Therefore, in our 
improved projection operator method, it is possible to examine the time evolution of conserved quantities which are 
composed of not only system operators but also environment operators. Any operator O{to) which commutes with 
the total Hamiltonian is conserved, for any order of expansion, because LOitg) = in Eqs. (|l3|) and (30). This is a 
reasonable result, because conserved quantities should be time independent. 

Until now, we have used the projection operator defined in Eq. (p^). However, to derive Eq. (pO|), only the condition 
( p^ ) is needed. Therefore, we can use a more general projection operator to derive Eq. (pO[) in the case that the 
condition ([19h is met. 



III. THE HEISENBERG EQUATION FOR A TIME-DEPENDENT HAMILTONIAN 

In this section, we consider the case of a Hamiltonian with explicit time dependence. We distinguish the time 
dependence of operators from explicit time dependence. To do this, the time dependent Hamiltonian H{t) is written 
with two time arguments: 

H{t)=H{t,t), (32) 

where the first argument represents the explicit time dependence, and the second that of operators. 
The time development of any operator O which has no explicit time dependence can be expressed as 

i f dsL{s,t(i) 

0{t) = e> O{to), (33) 

ij drL{r.tQ) 

where tg is the time at which we prepare an initial state. The operator e—y ° is defined as follows: 



I C drL{r,t„) ft fti .t,.-i 

^° =l + Y,i" dti dt2--- dtnL{tnM)L{tn-l, 

^ — 1 ^ to J to J tf) 



iJ dTL(T^to) 

Similarly, e^'" is defined as 



{t2,to) ■ ■ ■ L{tn, to). 



Therefore, the time evolution of the operator included in the Hamiltonian is expressed as 



(34) 



(35) 



i f dsL(s,tQ) 

H{t,t) = ed*° H{t,tQ). (36) 

By using this relation, it becomes clear that the operator 0{t) given in Eq. ( ^3|) satisfies the Heisenberg equation 

-0{t)^eJ"> iL{tM)0{to) 

i f dsL{s,tQ) 

= eJ'° i[H{t,to),0{to)] 

= i[H{t,t),0{t)]. (37) 

Note the following properties: 
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d iVdTL(T.to) , i r dTL{T,ta) 

as 



d i V dTL{T,to) i f* dTL{T,to) , 

— =e^= -iL(s,io , 

as 

(/if* drL{r.to) —i f* drL{r.tQ) d —i f * dTL(T,to) i dTL(T,to) 



(39) 



dt ds ^ ' 

From Eq. (|40|), we can derive the identities 

i r dr L(r^tQ) —i f dr L(T^tQ) —i f dr L{r.tf)) i f dr L{r.tQ) 

eJ' e^-'" ^e^-*' ' eJ' ' =1. (41) 

The following equations are derived similarly to those derived in the previous section. Corresponding to Eqs. (|^) 
and (^), we obtain 

ff if dsL^s^to) i f dsL{s.tQ) i f dsL{s,to) 

— eJ'" P = eJ'" PiL{t,to)P + eJ'° QiL{t,to)P, 

fl if dsL{s,to) i f dsL{s.to) i f dsL{s.to) 

-eJ'° Q = eJ'" PiL{t,to)Q + eJ'° QiL{t,to)Q. 

Equation ( ^3| ) gives us 

i dsL{s,to) i r dsL(s,to)Q f* i T dTL(TAo) i f* dTL(T,to)Q 

eJ'" Q = QeJ'° + / dse4'° PiL{s,to)QeJ' 

J to 

i f dsL(s,ta) /"* f dTL{T,ta) i f dTL(T,to) 

+eJ'° {P + Q) / dse^''" eJ'° 

Jta 

xPiL{s,to)QeJ^ 

i f dsL{s,tQ)Q i f dsL{s,to) 1 



where 



ft -i r drL{T.to) i r dTL(T.to) i f dTL{T.to)Q 

Sox(i,io)= / dse^ '° ed " PiL{s,tQ)QeJ^ . 

J to 



Therefore, we can obtain the following equation: 



(42) 
(43) 



(44) 



(45) 



W if dsL(s,to) 

0{t)=ed'° PiL{t,to)0{to) 

at 

i I dsL{s ,to) 1 

+e-:'° PSe,(i,to)- ———tL{t,h)0{to) 

t — ^0x1,1, to) 

i p dsL(s,to)Q I 

+geJ'° _____,L(^,i„)o(to). (46) 

This equation is equivalent to the Heisenberg equation, and corresponds to Eq. (13) in §2. When we ignore the 
explicit time dependence of the Hamiltonian, this equation agrees with Eq. (p^. 

Now we define the projection operator. As in the previous section, we consider the case in which the total Hamilto- 
nian can be divided into three parts, the system (S), the environment (E), and the interaction (I) between the system 
and the environment: 
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H{t,t)=Hs{t,t) + He {t,t) + Hi {t,t). (47) 

The self- interaction of the system and/or the environment can also be included in Hi(t,t). Furthermore, we assume 
the same form of the initial density matrix, given at the initial time to, as in Eq. (|l5|). We then define the projection 
operator as 

PO = Tte[peO] = {0)e (48) 
for any operator O. From the nature of the projection operator, we obtain the relations 

PLs{t, to) = Ls{t, to)P, QLsit, to) = Ls{t, to)Q, (49) 

LE{t, to)P = 0, LE{t, to)Q = LE{t, to), (50) 
where La{t, to)0 ~ [Ha{t, to), O] for a = S,E, I. Using these properties of the projection operator, we have 

QLo{t,to)Q = QLo{t,to), (51) 

i f (isLn(s,io) i [ dsQLQ{s,tQ)Q i f rfsLo(s,£o) 

QeJ^l) Q^QeJjt) , (52) 

where Lo{t, to) = Ls{t, to) + LE{t, to). 
Next, T,cx{t,to) is expressed as 

/•* -i P dTL{T,to) i P dTL{T,to) i r' dTL{T,ta)Q 

E^^(i,to)= / dse^ '° e^" PiL{s,to)QeA' 

J to 

— i f dsL{s,to) i f dsL{s.tQ)Q 

= Q~e^ ■'"> QeJ'" 
= Q- Uo\t,to)C,At,to)QV,^{t,to)Uoit,to) 
= Q- Uo\t,to)QUo{t,to) - UoHt,to){Ce^it,to) - l)QUo{t,to) 
-Uo\t,to)Q{V,^{t,to) - l)Uo{t,to) 

-U^\t,to){C,^{t,to) - l)Q(Pex(t,io) - l)Uo{t,to). (53) 
Here, we have introduced the following functions: 

Cox(i,io) = l + ^(-«)" / dti dt2--- dtn 

xLf^{ti^t())Lf^{t2,to) ■ • • Lf^{tn, to), 

V,^{t, to) = l + ^i" dti dt2--- dtnQLf °"(t„, to) 

xQL? '''{tn-i,to)---QLf '''{t,,to). 



where 



(54) 



(55) 



i f dsLQ{s,to) 

Uo{t,to) ^ e4'° , (56) 

Lrit,to) = Uoit,to)Li{t,to)Uo\t,to), (57) 
if "="(^,^0) = Uoit,to)Li{t,to)QUo\t,to). (58) 
By using the above functions and mathematical induction, we obtain 

i I dsL{s,tQ) 

= eJ'° PiL{t,to)0{to) 

i r* dsL{s,to) ]^ 

-e^'" PUo\t, to)C,^{t, to)Q— —————— C/o(i, to)iL{t, to)0{to) 

-L + [<~-cx[t, to) — ijC^ 

i r* dsL{s,to)Q I 

+Qe4'° ————iL{t,to)0{to). (59) 
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This is the final expression of the Heisenberg equation of motion with time-dependent Hamiltonian. Perturbative 
expansion up to first order in PEox(^, ^o)/(l ~ "^cxit, to)) leads to 

= ei '''''^''*°V{/-i(i,io)PC/o(t,to)iL(i,io)0(to) 

+ei'^'^'''*"V[/o-i(t,to)F / dsUo{s,to)iLi{s,ta)Uo\s,to)QUoit,to)iL{t,to)0{to) 



J to 

i r' dsL{s.to)Q 1 

+QeJ'« _____^L(i,io)0(io). (60) 

Here, we do not expand the third term on the r.h.s. of Eq. (|60|). This equation is the TCL equation in the case 
of a time-dependent Hamiltonian, and it becomes the same equation as Eq. (|3^) when we ignore the explicit time 
dependence in the Hamiltonian. The TC equation can be obtained as in Appendix C. 

IV. THE DAMPED HARMONIC OSCILLATOR MODEL 

In this section, we apply our improved projection operator method to the exactly solvable model called the damped 
harmonic oscillator model, [p^ which has no explicit time dependence. The Hamiltonian is 

H = Y, {MA-lAk - bIBu) + iTk{AlBl - AkBk)] . (61) 
fe 

The k dependence of is such that Tk is finite. We impose the following commutation relation: 

fe 

[A„ A]] = [B„ S]] = <5,,,, [A,,B,] = 0. (62) 
The time development of the operator j4^(i) is exactly solved as 

4(i) = e^'''* (cosh(rfct)4(0) + sinh(rfeOSfe(0)) . (63) 
This is a solution of the following equation of motion: 

|(4(t)) = in,{Al{t)) + rfctanh(rfct)(4(t)) + ^'''"-'-^k ^^^iVut) ^^"^"^^^ 

(64) 

Here, ( ) represents the expectation value taken with respect to an arbitrary density matrix p that satisfies the 
condition (p^. 

Now, we apply our method to this model, and compare the result with the exact equation (|64[). We consider the 
degree of freedom of the system to be A^ and that of the environment to be B^. The Hamiltonian is divided into 
three parts: 

Hs=Y,^kAlAk, (65) 

k 

He = -J^nkBlBk, (66) 

fe 

Hi=J2 *rfe(44 - ^kBk). (67) 

fe 

The projection operator is defined in Eq. (p^). When we substitute to = Q and 0(0) = A^ into Eq. (|3l|), we have 

|(4(t)) = ^n,{Aiit)) + Tlt{Alit)) + e^^-'r,{B,iO)). (68) 
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This equation is identical with Eq. ( pi] ) when we expand ta,nh{Tkt) and l/cosh(rkt) to lowest order in F,. Further- 
more, higher-order contributions are easily calculable. By using Eq. ( |D2| ) in Appendix D instead of Eq. (^), we can 
obtain the equation of motion with a higher-order term: 

|(4(t)) = zn^iAlit)} + rlt{Alit)) + e^^^'T, (^1 - ^rlf'^ (i?,(0)). (69) 

This equation also coincides with Eq. (^) with the appropriate pcrturbative expansion of F,. 

Now, we apply this model to the previous projection operator method. The approximate equation in that method 
which corresponds to Eq. ( |3l| ) is 

—0(t) = e^^(*-*")PiLO(to) 
dt 

J to 

+Qe'^Q(^~*o)ilO{ta). (70) 
Here, we use the same initial density matrix and projection operator. Now, the equation of motion for A^, is 

|(4(i)) = *0fc(4(t)) + F2t(AU0) + r.(s^(0)). (71) 

It seems that the factor e*^*"* in the third term on the r.h.s. of Eq. (^) cannot be reproduced. However, Eq. (|7^) gives 
consistent results only when we use a restricted density matrix that satisfies the condition PLe — 0. This condition 
implies 

PLsBkiO) = (Bfc(O)) = 0. (72) 

Therefore, the above equation is consistent with the exact result. From this argument, we can see that the previous 
projection operator method is valid only for a restricted initial state. 



V. CONCLUSIONS AND DISCUSSION 



We have given a new expansion of the Heisenberg equation of motion with a projection operator in the two cases 
of Hamiltonian both with and without explicit time dependence. In our method, one can prepare more general initial 
states which are forbidden in the previously applied projection operator method. Until now, the projection operator 
has been chosen as Eq. (|l6|) . However, if the condition QLqQ = QLq (in other words, PL^P = LqP) is satisfied, one 
can use Eqs. (|o|) and (|59|) for any kind of projection operators. 

Recently, a similar expansion method with a projection operator was proposed by Uchiyama and Shibata. In their 
method, one can expand the Heisenberg equation without such a restriction on projection operators as QLqQ = QLq, 
which is needed in our method. However, they impose some conditions on the operator O(to) whose time evolution we 
want to calculate. The first condition is that the time evolution with the unperturbed Hamiltonian should be solved 
as 

e^Loit-to}o(to)^fit,to)Oito), (73) 

where f{t, to) is a c-number. The second condition is 

POito) = O(io)- (74) 

This condition implies that we cannot calculate the time evolution of environment operators, when we use the definition 
(|l6|). However, it is possible to formulate the method without the above conditions. This extended Uchiyama-Shibata 
projection operator method is explained in Appendix E. 
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APPENDIX A: THE DEFINITION OF THE OPERATORS C AND V 



We define the operators C{t, to) and 'D{t, to) as 

C{t,to) = e'^o(*-*o)e-'^(*-*<'), (Al) 

^iQLQ{t-to) _ X>{t,to)e^''^'"°^^*~^°\ 

V{t,to) = e*'3'^^(*-*'')e"''^^°'^^*"*°). (A2) 
These operators satisfy the following differential equations: 

^C(t,to) = e'^°^*-*°\iLo - ii)e-'-^(*-*'') 
at 

= -iLi{t,to)C{t,to), (A3) 
^V{t, to) = e*'3i'Q(«-to)Q(j^ _ zLo)Qe-''3^oQ(t-to) 

= V{t,to)QiLf{t,to), (A4) 

where 

Li{t, to) = e'^°^*-''''>Lie-'^<'^*-'°\ (A5) 
L^(<, io) = e*^°(*-*'')L/Qe-'^°(*-*°). (A6) 
Prom these differential equations, we obtain 

C(t, to) = 1 + / '^^i / ^^^2 • • • / - ^0)^-/(^2 - to) 

X •••L7(i„-io), (A7) 

V{t,to) + dti dt2--- dtnQLf{tn-to)QLf{tn-l-to) 

n=l "^*o ''to "'to 

x---QLf{ti-to). (A8) 
Similarly, Cex(t, to) and X>ex(t, to) are defined as 

— i d5L(s,fo) 

e^-^'" =C/o-\t,to)Cex(t,to), 

— t r dsL{s,to) 

CUtM) = Uo{t,to)e^ " , (A9) 

i f dsQL{s,to)Q i f d.sQLo{s,to)Q 

eJ*" =2?ex(t,to)eJ"' 

i f dsQL{s,to)Q —i f dsQLo{s,to)Q 

V^^{t,to) = eJ'° e^-"*" . (AlO) 

Differential equations can be constructed again, and are solved to yield 

Cex(t, to) = 1 + / dti f ' dt2-- - [ " ' dtnLT{ti,to)LT{t2, to) 

x---Lr(t„,to), (All) 

2?ex(t,to) = l + / a!ti / dt2--- dt„QL^°"(t„,to)QL?'"(t„_i,to) 

x---Q£^^''(ti,to), (A12) 



where 



Lr(t,to) = Uo{t,to)Li{t,to)U^\t,to), (A13) 
Lf '"(t,to) = C/o(t,to)L/(t,to)0[/o-'(t,to). (A14) 
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APPENDIX B: THE TRANSFORMATION OF THE OPERATOR S(r,To) 



The operator S(t, to) in §2 can be expressed as 

Jtn 



d 



to 



ds 



f p _ ^-iL{t-to) p^iL(t-to)Q _ r ^g^rf,-iL{t-s) ^iL(t-s)Q\\ 



'to 

— — {P — Q-iHt-to) p^iL{t-to)Q g-si(t-to)gji(t-to)<3| 

= Q — g-iL(t-to) QpiL(t-to)Q 

= Q- e-'^'>^'-^'>^C{t,to)QV{t,to)e''^^'''^^*~*°^. 
Similarly, Sex(i, io) in §3 is given by 

ft -i f* dTL{T,to) i r dTL{T,to) I f dTL(T,to)Q 

Sex(i,to)=/ dse^ '° eJ*° PiL{s,to)QeJ^ 

J to 



i 



* ~* I! '''^-^('^'*o) * // <lTL(T,to) d i r* dTL{T,to)Q 

dse^ ° P—e^" 
to ds 

—i r* dsL{s,to) i f dsL{s.,to) i f* dsL{s,to)Q 

*» {gj'o P-PeJ*° 

ft i(' dTL(T,to) ^ ^ i r dTL{T,to)Q, 

/ dseJ'° iL{s,to)PeJ' } 

J to 

—i r* dsL{syto) r ^ r* dsL(s,to) i f dsL{s,to)Q 

*° JgJ*o P-PeJ*° 

d ( <*«^(«'*o) i f ' dsL{s,to)Q\ 1 

-i f dsL{s.,to) (■ i f dsL{s,to) i dsL{s,to)Q 



i r* dsL{s,to) i f dsL{s,to)C 

''*0 _l_ X, *o 



} 



—i f dsL{s,to) i f dsL{s,to) —i f* dsL{s,to) i f dsL{s.,to)Q 

—i r* dsL{s,tQ) i f dsL[s,to)Q 

= Q-e^-''° QeJ*° 

i r' dsQLo{s,to)Q 

= Q- Uo\t,to)Ce^{t,to)QV^^{t,to)eJ*° 



To allow for simultaneous discussion, we introduce the following notation: 

C{t) 



p-»Lo(t-to)c(t, to)e^-f^°(*-*<') for the case of § 2 

UQ^{t,to)Ccx{t,to)Uo{t,to) for the case of § 3, 

g-»Lo(t-to)p(i^ ^^)g<Lo(t-to) for the case of § 2 

[/o"^(i,to)2?ex(i,to)C^o(i,io) for the case of § 3, 

g-iLo{t-to)Q^iLo{t-to) for the case of § 2 

UQ^{t,to)QUo{t,to) for the case of § 3. 
Using the mathematical induction, we confirm the following relation: 



m = 

Q{t) = 
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pi:{t,to){Qi:{t,to)r 

- [(-l)"-ip{QW(C(t) - i)}"g(t) + (-i)"-ip{(c(t) - i)Q(t)}"+i] 

n-l 

-I 



+pJ2{-iy{Q{tmt) - i)yQ{t){v{t) - i){Qi:{tMr-'- 

1=0 
n-1 

n 

-pY,{-iy{Q{t)ic{t) - i)yQ{t){v{t) - i)(QE(t,io)r-' 

n 

-pJ2{~iy{m - i)Q{t)y+\v{t) i){Qj:{tMr-\ m 



1=0 



where n is integer and n>l. The second and third terms in PI](t, to)(QS(t, to))" and the fourth and fifth terms in 
PE(t, to)(QS(t, to))"~^ cancel. The fourth and fifth terms in PE(t, to)((3S(t, to))" and the second and third terms 
in PE(t, to)(QE(t, to))""*"^ also cancel. Therefore, only the first term survives. As a result, all the terms including 
!?(<) disappear. Finally, noting the relation E(t,to) = E(t,to)Q, we find that PE(t, to) j^rj^^n^ '^^^ expressed as 
follows: 

PS(t,to)- — - = PE(t,^o)- 



l-E(t,^o) ^ ' '^^l-QE(t,to) 

oc 

= PS(t,to)E(QE(t,to))" 

oc 

= -PY.[{-Q{tmt) - i)}"Q(t) - {-(c(t) - i)g(t)}"+i] 

n=0 

= -PQ(t) ,,,,-P. (CW-l)QW 



1 + (C(t) - l)Q(t) 1 + (C(t) - l)Q(t) 



APPENDIX C: THE DERIVATION OF THE TIME-CONVOLUTION EQUATION 

Here, we derive the equation with a time-convolution integral. In this appendix, we consider the case of a Hamil- 
tonian without explicit time dependence. The discussion, however, can also be applied to the case of a Hamiltonian 
with explicit time dependence. We substitute Eq. (^, (instead of ([III)) into Eq. (|^), and we obtain 

-rO{t) = e'-^(*-*«)PiLO(to) + / dse^^^"-^"^ PiLQe'^^^*-'hLO{ta) 



dt 



+Qe*^'3(*^*°)iLO(to) 



e'^^*-*o) PiLO{to) + / dse'^^*-'^ PiLQe'^^^'-'"hLO{to) 

J to 

+Qe'LQit-to)-^Q(^^^y (CI) 



This equation is equivalent to the Heisenberg equation. For the existence of the operator e*^(*~'*\ the second term 
on the r.h.s. of the equation has a time-convolution integral. Therefore, this equation is a called a time-convolution 
(TC) equation. This can be rewritten by using X'(t,to) as 

—0{t) ^ e'^^'-'"'^ PiLO{tn) + / dse'^^'-''>PiLQV{s,tn)e''^^°'^^''~'°hLO{to) 



dt 

+Qe'^Q(t-*o)iioito). (C2) 
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By expanding V{t,tQ) to lowest order, we have|^ 

''-O(t) = e'^'-'-'"^ PiLO{to) + I dse'^^'-'^PiLQe''^^°'^''''-'°hLO{to) 



dt 

+Qe'^<3(*-*°)jLO(to). (C3) 

At this level of expansion, the difference between the equation in the improved projection operator method and that 
in the previous projection operator method cannot be observed, because this difference is included in the higher-order 
term in 'D{t, to). 

In the case of the TC equation, the expansion ( |C2| ) is always valid, because no restriction is imposed on the 
projection operator. This is different from the case of the TCL equation, for which the condition QLqQ = QLq must 
be satisfied. 

APPENDIX D: THE DERIVATION OF HIGHER-ORDER TERMS 

We expand C{t,to)/{l + {C{t,to) — 1)} to second order in the interaction: 

1 



C{t,to)Q 



l + (C(i,to)-l) 

C{t, to)Q{l - (C(t, to) - 1)Q + iC{t, to) - l)Q{C{t, to) - 1)Q + • • ■} 



t rt rt 



dtiLi{ti-to) + {-t)^ dti dt2Li{ti-to)Li{t2~to)} 



to J t.Q J to 

x{l + i dtiLi{ti-to)Q - i-i)^ dti dt2Lj{ti^to)Lj{t2-to)Q 

J to J to to 

+ {-i)^ f dtiLi{ti-to)Q f dt2Li{t2-to)Q} 

pt [■t ptl 

Q-P dtiiLi{ti-to)Q + P dti dt2iLi{ti~to)PiLi{t2-to)Q 

Jto J to J t.Q 

-P f dti f dt2iLi{t2-to)QiLi{ti-to)Q. (Dl) 

Jto Jto 



In substituting this into Eq. (pQ|), we have 



_^g«L(t-to)pg-»Lo(t-io)p f dsiLi{s - to)Qe'^°^*-*°hLO{to) 



Jto 

_^iLit-to)p^-iLoit-to)p I ^^f' dt2iLi{ti ~ to)PiLi{t2 - to)Qe'^°^*-'°\LOito) 

Jto Jto 

^gzL(t-to)pg-zLo(t-to)p f /" ' dt2iLi{t2 - to)QiLi{ti - to)Qe'^°^'-'°\LO{to) 

Jto Jto 

+ge'^Q(*-'») ^_^^^^^^^^ »LO(to). (D2) 
The third and fourth lines on the r.h.s. of this equation are higher-order contributions. 



^It may seem difficult to solve e'"'^'^"'^^" ^°''iLO{to), because of the operator Q. However, it can be solved using the following 
differential equation: ^e*'3^o'3(=-'o)iLO(fo) = e''^^°'^''^-'°\QLoQiLO{to). 
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APPENDIX E: THE EXTENDED UCHIYAMA-SHIBATA PROJECTION OPERATOR METHOD 



Here we explain the extended Uchiyama-Shibata (U-S) projection operator method. The time evolution of an 
operator can be written as 



where 



i f dsL{s,to) 

0{t) = pj'« Oih) 

= U-{t,to){Uo\t,to))^0{to), 

i f dsLi{s,tQ) 



(El) 



(E2) 



Here Li{t, to) is the same as that defined in §2. Now, we expand all the fields with creation and annihilation operators 
at to, and prepare the initial state at to- The time evolution is determined by the following equation: 



|o(t) = |;7_(<,io)-({/crHt,M)^0(to) + ?7-(i,to)|({/o-'(i,M)^0(io). 



(E3) 



Let us now discuss the time evolution of lJ-{t,to). We introduce a generic projection operator which has the 
following properties: 



p2 = p_ 

Q = i-P, 



(E4) 
(E5) 
(E6) 



These are the same as those introduced in §2. By using the projection operators, we obtain the following differential 
equation: 



dt 



U-{t,to) - U-it,tQ)PiLi{t,to) + U-^it,to)QiLi{t,to). 



We can construct the differential equation for U^{t,to)Q and solve it as we have done in §2: 



U-{t,to)Q Qu-{t,to) + / dsU-{s,to)PiLi{s,to)Qu^{t,s) 

J to 

= [Qu-{t,to) - U-{t,to)P{ez\t,to) - i)]e-{t,to), 



where 



f dsiLji 



s,to)C 



u-{t,to) = e- 
U+{t,to)^e. 
e_(t,to) = [l- / dsU+{t,s)PiLi{s,to)Qu^{t,s)]^^. 



- r dsiLi{s,tQ) 

JtQ 



(E7) 

(E8) 
(E9) 

(ElO) 
(Ell) 
(E12) 



When we substitute Eq. (EJ) into (E7), we can obtain the TC equation. On the other hand, substituting Eq. (E9) 
into (E7), the TCL equation can be derived. Here we discuss only the TCL equation. We obtain 



dt 



U-{t,to) = U-{t,to)PiLi{t) - U-{t,to)P{l - Q-{t))iLi{t) + Qu-{t,to)&-{t)iLi{t). 



(E13) 



When we substitute this result into Eq. (E3), we obtain 
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dt 



+Qu-{t,h)e-{t,to)tLi{t,to){U^\t,to)yO{tn). (E14) 



This equation is equivalent to the Heisenberg equation of motion and corresponds to Eq. (13). Up to second order, 
we have 

-0{t)^e4'° iLoO{to) + e^'" {Uoit,h))^ PtLi{t,to){U^\t,to)y Oito) 

+eJ'°''''^'^'''"\uo{t,to)yP [ dsiLj{s,to)QiLi{t){Uo\t,to)yO{to) 

Jto 

+Qu-{t,t„)&^{t,to)tLi{t,to){Uo\t,to))^0{h). (E15) 



Here, we do not expand the third term on the r.h.s. of Eq. ( E15 ). 

The derivation given here is not the exactly same as that proposed by Uchiyama and Shibata. We have modified 
the following two points. In the original paper of Uchiyama and Shibata, they imposed the condition 

iU^\t,to)yO(to) = f{t,to)0{to), (E16) 



where f{t,ta) is a c-number function, and O{to) is the operator seen on the r.h.s. of Eq. ( |E14| ). However, we have 
not used this condition in the derivation given here. Furthermore, they imposed the condition 

POito) = Oito). (E17) 

It is easily seen that no condition needs to be imposed to derive the extended equation in the extended U-S method, 
while in the method given in this paper the condition QLqQ = QLq is required. However, there is a merit to our 
method: All quantities which commute with the total Hamiltonian are conserved in the time evolution regardless of 
the order of the expansion of the equation. It is not clear that this property is satisfied in the extended U-S method. 

This can be shown concretely. We consider the following Hamiltonian: 

H = Uaa^a + Ubb^b + g{aH'^^ + a)H'^). (E18) 

The operator a is the degree of freedom of the system, and the operator b is that of the environment. The total 
Hamiltonian can be divided into the following three parts: 

Hs = cJafl^a, (E19) 
He = LUbb^'b, (E20) 
Hi = g{aH^^ + a^^^). (E21) 

We use the same projection operator as in Eq. (p^. Now, we calculate the time evolution of the total Hamiltonian 
H in the extended U-S method. We expand (1 — G)_(t)) to zeroth order in the interaction in Eq. (E14). This means 
that we ignore the third term on the r.h.s. of the equation. The equation of motion is expressed as 

— {H(t)) = (e'^HLoH + e'^'e~'^°*Pe'^°HLiH) 
dt 

= 2ig{^a - ujb){-{a\t)b^\t)) + {a'^\t)b\t)) 
+e2*"''*(a2(t)}(6t2) _ e-2»"''*(at2(t))(62))^ 

(E22) 

where we take t^ — 0. The total Hamiltonian has no time dependence originally. However, it is not obvious whether 
the r.h.s. of the equation becomes zero for arbitrarily states. 

On the other hand, the time dependence of the Hamiltonian calculated using our method vanishes: -^H{t) — 0. 
This can be easily shown. To begin, we substitute O — H into Eq. (|3l|). Then, we must calculate the commutation 
relation LO = LH first. This quantity clearly becomes zero. This is satisfied regardless of the order of the expansion 
and Hamiltonian considered. Similarly, other quantities that commute with the total Hamiltonian are also conserved. 
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